function f=Tmj_Normal(x,m,j,mu,sigma,a,b)

if (m==1)
    f=T1j_Normal(x,j,mu,sigma);
end

if (m==2)
    f=-(b+a)/(b-a)*T1j_Normal(x,j,mu,sigma) + 2/(b-a)*T1j_Normal(x,j+1,mu,sigma);
end

if (m>=3)
    f=4/(b-a)*Tmj_Normal(x,m-1,j+1,mu,sigma,a,b) - (2 + 4*a/(b-a))*Tmj_Normal(x,m-1,j,mu,sigma,a,b) - Tmj_Normal(x,m-2,j,mu,sigma,a,b);
end

    %
